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The homotopy theory of topological defects is a powerful tool for organizing and unifying many 
ideas across a broad range of physical systems. Recently, experimental progress has been made 
in controlling and measuring colloidal inclusions in liquid crystalline phases. The topological 
structure of these systems is quite rich but, at the same time, subtle. Motivated by experiment 
and the power of topological reasoning, we review and expound upon the classification of defects in 
uniaxial nematic liquid crystals. Particular attention is paid to the ambiguities that arise in these 
systems, which have no counterpart in the much-storied XY model or the Heisenberg ferromagnet. 



flux lines is essential for maintaining superconductivity 
(Nelson 1988^, and the theory of magnetic monopoles 
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allows them to be viewed as point particles (Coleman 



1988). In all of these cases, the fluctuations and configu- 



First identified through the beautiful textures of their 
defects, liquid crystalline materials may very well be 
the ideal proving grounds for exploring notions of bro- 
ken symmetries, associated Goldstone modes, phenom- 
ena akin to the Higgs mechanism, and, of course, topo- 
logical defects. Indeed, the subjects of topological de- 
fects and liquid crystals are so intertwined that it is dif- 
ficult to see how a thorough understanding of one could 
be garnered without knowledge of the other. Moreover, 
it has proven fruitful to treat the defects as fundamen- 
tal excitations; solitons in the sine-Gordon model can 



rations of the smooth background field around the defects 
can be replaced with an effective interaction between the 
defects. Thus the elasticity and statistical mechanics of 
these systems can be recast in terms of a discrete set 
of topological charges. Because of these successes and 
physicist's deeply-ingrained love for Gauss's law and the 
connection between flux and charge, we naturally at- 
tempt to interpret more complex systems in the same 
terms. Unfortunately, these two cases are exceptions - 
in this colloquium we will emphasize this important point 
and attempt to clarify if, when, and how it is appropriate 
to view defects as point charges. 

Nematics are a prototypical liquid crystalline material 
- though a sample will flow like a liquid, owing to the 
lack of positional order, the rod-shaped, typically organic 
or aromatic ring molecules are locally aligned, leading 
to orientational order, with widely-applied optical conse- 
quences. Defects in a nematic liquid crystal, as in any 
ordered medium, are places where the order changes dis- 
continuously and is thus ill-defined. They occur in the 
form of isolated points and lines, either system spanning 
with endpoints on the boundaries of the sample or closed 
up into loops. We will use as our example the three- 
dimensional uniaxial nematic, a system that turns out 
not to be as simple as it might seem. Not only is it, in 
some sense, the simplest counterexample to the "defects 
as charges" approach, but it is also now extremely rele- 
vant from the experimental point of view as it has been 
demonstrated that defects can be manipulated in nematic 
cells with great precision and variety allowing complex 



be treated as interacting fermions (|Coleman[[l975[), vor- hnks and knots to be tied ( |Kamien[ |2011[ |Tkalec, et al. 



tices in the XY model drive the Kosterlitz-Thouless phase 
transition (Kosterlitz and Thouless 1973), controlling 



2011). Though the history or path dependence of the 
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defect motion at once spoils the program of reducing the 
configuration space to discrete data, it simultaneously 
introduces a new set of topological degrees of freedom 
in the three-dimensional nematic. The path-dependence 
imposes a non-Abelian structure to the classification of 



states. In other liquid crystalline systems, line defects in 
the bi axial nematic are anticipated to entangle topolog- 
ically (IPoenaru and Toulousel 119771) and the interaction 



between disclinations and dislocations in smectics suffers 



similar path dependence ( Chen, et al. 2009^ Pocnaru and 



Toulouse 1977). Certain non-Abelian quantum systems 



have been proposed as a route to topological codes and 



computation (Kitaev 2003 Nayak, et al. 2008) 



Although they arise naturally during phase transitions, 
defects can also be induced deliberately by local exci- 



tation with a laser ( Smalyukh, et al. 2009 ) or through 



boundary conditions as with colloidal inclusions where 
the anchoring of molecules at the colloid surface induces 
the presence of defects in the surrounding liquid crys- 
tal. This is an especially rich technique allowing for the 
generation of a wide range of defects and for studying 
the interactions between them. For instance, the defect 
accompanying a colloid may take the form of an isolated 
point 



cling t 



Poulin, et al. 1997) or of a disclination loop encir- 



le particle in a 'Saturn ring' configuration (Teren- 



tjev 1995), while two or more colloids in close proximity 



can form a variety of 'entangled' structures where a single 
disclination loop wraps around both of them (j Guzman, et 
aLJ |2003||Ravnik, et aL[|2007||Ravnik and Zumern2009p . 



These last situations realize the interaction between line 
and point defects that is a central aspect of this collo- 
quium. 

In Section [n] we begin our tale by reviewing the classi- 
fication of topological defects by looking at point defects 
in two-dimensional nematics, and line defects in three- 



dimensional nematics. In Section III we continue to point 
defects in three-dimensional nematics. With all the char- 
acters introduced, we proceed to discuss our main point: 
that there are only two types of unlinked disclination 
loops. What we present is built on the classic review 



by Mermin on topological defects (Mermin 1979) and 



a definitive result by Janich (Janich 1987) which finds 
the topological classes of uniaxial nematic textures on S'^, 
the three-dimensional sphere, (equivalently, three dimen- 
sional space with fixed boundary conditions at infinity), 
with point and line defects. We will also discuss the bi- 
axial nematic phase and show that it affords additional 
physical insight into the analysis of uniaxial defects. 



II. CIRCULAR LOOPS AROUND DEFECTS: THE 
FUNDAMENTAL GROUP 

A common and distinctive experimental procedure for 
imaging the orientational order is to place the sample 
between crossed polarizers, yielding in nematic materials 
the characteristic Schlieren texture from which they de- 
rived their name. In FIG. [T] our eyes are drawn to the 
dark brushes and tend to follow them towards their inter- 
sections. The brushes typically meet at points in twos or 
fours; larger numbers are possible but they too are mul- 
tiples of two. The dark brushes correspond to regions 
where the orientation is parallel to either of the mutually 



perpendicular polarizer or analyzer directions and the 
points at which they meet are disclination defects. Why 
defect? Because there the orientation is ill-defined, as the 
brushes would tell us to assign two or more different ori- 
entations to that point. We can fix this by removing the 
defect point by poking a hole in the material. Physically, 
of course, we don't have a hole in our sample, just a place 
where we don't know the order. The problem is fixed by 
letting the magnitude of the order vanish at the origin. 
Then no hole is needed, but there is a point where there 
is no longer nematic order and the core is in the isotropic 
phase. To the reader acquainted with defects in super- 
fiuids and superconductors, this should sound familiar: 
in the center of an Abrikosov vortex in a superconductor 



(Abrikosov 1957) the superconducting order parameter 
vanishes and there is normal metal - a hole in the super- 
conductor (Chaikin and Lubensky 1995). In all these 



cases we can detect the existence of a defect without any 
detailed knowledge of the disordered core: Topological 
defects are characterized by the boundary conditions at 
the interface between the higher symmetry state and the 
lower symmetry state. In the Schlieren texture, these 
boundary conditions may be read off from the pattern of 
brushes. 

In this section we introduce the basic concepts in the 
study of topological defects in this simplest case and con- 
nect the observation of brushes to the presence of defects. 
We will cover what defects are, the use of measuring sur- 
faces to probe them, and finally how to combine them. 
We begin by recalling that the continuum description of 
uniaxial nematics is based upon the assignment of a lo- 
cal 'average molecular orientation' to every point in the 
sample, called the director field. This average orientation 
is that of a rod-like object, rather than an arrow, so that 
the nematic director is properly a line field, rather than 
a vector field. Moreover, it is only the orientation of the 
director that is important, not its magnitude, so that we 
may take the director to be a unit vector n, subject to 
the condition that n and — n are identified. 



A. Two-dimensional Director: RP 

First we consider nematic samples in which the direc- 
tor always lies in the xy-plane, a prototypical thin cell 
situation. Equivalently, we consider the projection of a 
three-dimensional director onto the xy-p\ane: 

p = A{x, y) [cos {e{x, y)) , sin {e{x, y))] , (1) 

where A{x, y) is the nonvanishing magnitude of the pro- 
jection and 9{x,y) is the angle it makes with the x-axis. 
Just as n and — n are identified, so too are p and — p or, 
equivalently, 6{x, y) and 9{x, y)-t-7r. Because the polarizer 
and analyzer directions are 7r/2 apart, each dark brush 
around a defect point tells us about a 7r/2 rotation of 
the director. The number of brushes counts the winding 
of the director field. The fact that we have defects with 
two brushes directly tells us that rotations by tt leave the 



phase invariant and is experimental proof of the symme- 
try of the nematic phase (Frank 1958). The angle that 



the director turns on a loop around a defect divided by 
the angle 27r of a full turn, is commonly referred to as 
the strength, or charge, of the disclination and is equal 
to ±1/4 times the number of brushes. Though the abso- 
lute value of the strength may easily be extracted from 
counting the number of brushes, the sign of the winding 
is not yet determined. It is a pleasant exercise in visu- 
alization to show that upon rotating the polarizers the 
brushes will counter-rotate when the defect has negative 
winding and will co-rotate when the winding is positive. 
As defects are points of discontinuity of the director 
field, we aim for a classification up to continuous defor- 
mations of their surroundings. When two paths can be 
continuously distorted into each other, they are said to 
be (freely) homotopic. It is useful to think of these de- 
formations as time evolution because it not only allows 
us to visualize the distortion but it also reminds us that 
under continuous dynamics, configurations remain in the 
same homotopy class. Above we began to think about 
a defect in terms of a measuring loop around it - this 
key idea is what allows us to probe the properties of a 
singularity in our texture without leaving the safety of 
the well-behaved director field. Let's consider the possi- 
bilities for the director on a small round measuring circle 
in the system. Here's an explicit family of possible mea- 
surements on this circle p/A = c,„ which exhibit any 
even number of brushes: 



= \cos{m(j)), sin(m0)] , 



(2) 



with m G |Z a half-integer and the polar angle around 
our circle. These textures are all distinct in that we can- 
not continuously deform c„ into c„j for n ^ m. Thus 
within this family, the half-integer m classifies the pos- 
sibilities for the director field. Of course the textures 
described by ^ are special choices, or simplifications, 
and the director is rarely of this exact form even very 
close to the defect line. However, any director field that 
we might have measured on our circle can be homotoped 
into one of these! The utility of (l2|, therefore, is that 
these fields capture the possible windings of a general 
director field, and thereby form a set of representatives 
for the behavior outside point defects. Most placements 
of our measuring circle will not surround a defect and 
hence the director field can be deformed to Cq, which is 
a constant. On the other hand, if the measuring circle 
surrounds a point where two brushes meet, it should be 
homotopic to Ci/2 or c__i/2. 

Though reasonable sounding, the previously stated 
connection between brush counting and winding num- 
ber relies on a hidden assumption! Suppose p = 
[cos /((/)), sin /(0)] on our measuring circle and / took 
the form 



/l-^) 
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27r 
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(3) 
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FIG. 1 Cartoon of the director field and the resulting brush 
pattern corresponding to the homotopy in equation (4). The 
radial coordinate t parametrizes the sequence of maps ft{4i)- 
Near the defect center, six brushes are visible so that simply 
counting would give a winding of |. However, away from the 
defect core, the texture smoothly deforms to a texture that 
may be seen to have a winding of |. 



where we first go from to tt, back to and then back 
to TT. Now / takes values in the interval [0, tt] with and 
TT identified due to the p — > — p symmetry. Were we to 
count the number of times we landed at 7r/2 for instance, 
we would find 3, as we would for any generic point of 
[0, tt]. Under crossed polarizers we would see six brushes. 
But we can also smoothly distort the map /(</>) so that 
there are only two brushes through the sequence of maps 
ftW,te [0,1]: 



ft{4>) = 




2)tt 



0< 

3-2t 
(4-2t)7r 
3-2t 



27r 
3-2t' 



) < 

<(f>< 

<6 



(4-2t)7i- 
3-2t ' 

< 27r. 



(4) 

We see that /o(0) = /(0) while /i(0) generates the 
charge 1/2 map Ci/2 as ranges from to 27r. So have 
we shown that charge 3/2 is the same as charge 1/2? 
Certainly not! We have learned, rather, that we must 
be more careful about counting brushes. We assumed in 
our previous rule of counting brushes to determine angles 
that the angle function was always either increasing or 
decreasing on our loop. In general, when counting the 
number of times we pass 7r/2 we must keep track of the 
way we pass it - is the angle increasing or decreasing? 
We thus look at the net rotation of /((/)) by calculating 



the winding number via the following integral 






(5) 



for all t. In general, the winding number is precisely the 
charge of the defect. In actual liquid crystal systems, 
energetics strongly disfavors extra brushes because these 
would correspond to regions of large splay or bend in 
a nematic. Thus, the naive counting of brushes works 
outside of exceptional circumstances. Though we have 
framed the discussion above in terms of a particular mea- 
suring circle, the classification by winding numbers holds 
for any measuring loop disjoint from the defects - just 
let (f> be any coordinate which goes from to 27r as we 
traverse the loop once. In particular, the texture on any 
loop in the sample which lassos just one of the defects 
is homotopic to any other; this can also be seen by the 
fact that the number of brushes meeting at the defect is 
conserved. Similarly, the number of brushes meeting at a 
point in a Schlieren texture remains unchanged in time as 
the texture evolves except when two defects combine and 
annihilate each other. The winding number m is a key 
example of an essential idea in topology: turning geomet- 
ric information into counting. By doing so we get robust 
measures since, in particular, integers (or half- integers), 
being discretely valued, cannot change continuously. 



B. Making Precise Measurements 

Let's start putting our above classification into a 
broader context. We have mentioned several different 
topological spaces which we should keep straight. First, 
we are dealing with uniaxial nematic textures, so there is 
a space associated with the possible ground states of the 
director field, the ground state manifold (GSM)^. Here, 
the angle of the director T lives in the interval [0, tt] with 
and TT identified. This space is known as the real pro- 
jective line, RP"'^. Like any interval of the real line with 
endpoints identified, it has the topology of a circle. Next, 
we have the space M making up the sample volume. For 
our Schlieren textures, M may be taken to be a suitable 
region of the plane. Inside M there is a subspace, the set 
of defects E where the order is not defined. In the last sec- 
tion, E was the set of disclination points. We obviously 
have the director field, an MP -valued field on M away 
from the defect set, n : Af \ E -^- GSM. In general, it is 



^ Recall that the ground states of a nematic are just uniform tex- 
tures with the director pointing in some direction, so that in 
a general texture the local orientation can always be identified 
with one of these ground states. Thus the changing orientation 
of the director in a general texture may be thought of as motion 
on the ground state manifold. There is not an accepted term 
for this space in the literature, where it has also been called the 
manifold of internal states and the order parameter space. 



daunting to contemplate the texture on the entire sam- 
ple especially if there are many defects. We also would 
like to understand to what extent n may be understood 
as arising from its behavior near the defect set - might 
there be a way to cut M \ E into more manageable pieces 
around its defects, classify those, and then glue them 
back together? These considerations motivate the intro- 
duction of auxiliary spaces: measuring surfaces of fixed 
topology contained in the sample. We study the behavior 
of the director field on these as we vary their placement 
relative to the defect set or as we vary n. Indeed, all 
the topological "charges" in this paper are defined with 
respect to a choice of measuring surface. Saying that a 
defect carries a charge is really shorthand for saying that 
a small measuring surface which surrounds only that de- 
fect measures such a charge. If such a measuring surface 
is not present, then the "charge of a defect" is ill-defined. 
Of course, these ambiguities are nothing new and have 
been discussed before (Mermin 1979). However, with 



the interest in nematics and cholesterics with embedded 



colloids ( Copar and Zumer 2011 Lintuvuori ei aL 2010 



Musevic, et al. 2006), these mathematical subtleties are 



no longer just about precision or axiomatic rigor - they 
are absolutely necessary for the proper interpretation of 
data. Above we measured the director on loops (which 
are topologically circles, denoted by §^) and observed 
that the winding number m was constant during defor- 
mations of the loop or the texture as long as no defects 
pass through the loop, which essentially solved our clas- 
sification problem for single defects in Schlieren textures. 
In more mathematical jargon, we saw that the set of maps 
from S^ to RP^ up to homotopy (denoted [§\MP^]) was 
equivalent to the set of half-integers |Z. 

How can we go from single defects back to the full tex- 
ture? Let's draw a picture of a two dimensional plane 
with punctures at the defects and small measuring loops 
about each puncture. Given the winding numbers at the 
loops, do we know the classification of the full texture? 
This picture is actually no different from one we might 
draw for the following situation: Consider a hypothetical 
Schlieren texture that contains several nearby strength 
±1/2 defects. If we take a measurement around a circuit 
that surrounds all the defects, what do we get? Before 
we say anything more about the local-to-global question, 
let us try to find an answer to this natural question: how 
can the theory capture the intuitive notion that defects 
can combine or split? We almost have it: we know that 
measuring circuits around defects are in correspondence 
with half-integers, and it happens to be true that any two 
half-integers can be added or subtracted. There is just a 
tiny gap between question and answer now. For one, we 
need to define a way to go from two measuring circuits 
to one. The natural way to add loops is to concatenate 
them by forming a longer loop by running over the first 
and then running over the other. But in order to do this, 
the loops have to start and end at the same point in sam- 
ple or they cannot be connected. By the same token, the 
two textures traced out on the loops must start and end 



at the same point in RP ! Given two arbitrary measur- 
ing loops, there's no reason for the textures on them to 
agree at some point. On the other hand, our figure shows 
that we can draw a big circle surrounding any pair. Is 
there a way to go from two circles to one without having 
to make arbitrary choices in the region between the cir- 
cles? Yes! Since we are interested only in properties that 
are preserved under continuous deformations, we can de- 
form the texture around one defect so that the value in 
KP agrees at the common base point in the sample. To 
visualize this, consider a small disc around the defect. 
The homotopy proceeds radially outward by deforming 
the original texture at the center of the disc to the new 
texture at its outer edge. 

This solution is inspired by the following geometrical 
fact: the space between a set of non- intersecting circles in 
the plane can be continuously deformed to a set of circles 
joined at a single point, called a wedge sum, or bouquet 
of circles. Once we have the texture on the bouquet, the 
concatenation of the loops is well-defined, and we thus 
have a definition of addition in our system. In the case 
here, we may form addition by measuring around the loop 
which hugs the "outside" of the bouquet. By restricting 
to sets of loops which pass through the same base point, 
our set of equivalence classes of maps from §^ to the GSM 
is enriched with the additional algebraic structure of a 
group. Note that based homotopy is absolutely necessary 
here - continuous distortions of paths that preserve the 
base point. 

A quick aside on group theory is in order. Recall that 
a group requires a way of adding its elements: we just 
saw that this is the concatenation of loops. In terms of 
brush counting, or the winding number, we need not be 
concerned with the details of the rate at which we tra- 
verse the GSM, or whether we pass a brush in the first 
or last half of the trip. We need only concern ourselves 
with the order in which we put the two paths together. 
This addition of classes of paths defines the group addi- 
tion which may not be commutative but happens to be 
in this particular case. A group requires an identity: in 
a uniform nematic texture the map from a loop in the 
sample to the GSM is always constant, so this provides 
the identity under the group addition. Finally, each map 
has an inverse: because we can go around the GSM in 
the opposite direction by backtracking our precise path, 
we also have an inverse. Together these properties ensure 
that we have a group, known as the fundamental group, 
TTi (GSM) - the set of based maps from a measuring cir- 
cuit §^ in the sample to the GSM which are equivalent up 
to based homotopy. And as the reader probably has al- 
ready guessed, this group for the case where GSM — MP^ 
is precisely the (half-)integers. A loop that goes around 
only one defect in our hypothetical texture will intersect 
two brushes, giving a winding of ±1/2. A larger loop 
encircling both will yield a count of four if the defects 
have the same sign, or zero if they have opposite signs. 



Though the choice of the base point is arbitrary, its 
constancy is essential and leads to many of the interest- 
ing phenomena we will discuss. The need for a base point 
is precisely why standard time zones were developed for 
trains. It is just a matter of having a single clock at one 
location (i.e. one measuring circuit) to determine the 
elapsed time at that spot. Two clocks can be used to 
measure elapsed time as long as they are synchronized - 
that is why we must have a base point. How are general 
maps from S^ to the GSM related to 7ri(GSM)? In the 
case of Schlieren textures, we can relax: [S^, MP^] is sim- 
ply the underlying set of 7ri(MP^). Even here, however, 
we must remember that if we want a topological charge 
for loops around defects that satisfies addition, we better 
pin our loops on a base point! 

Let's conclude this section by sketching how what we 
have done tells us how to go from local measurements 
around defects to the full texture. The problem of clas- 
sifying the defects on the full texture if we fix a base 
point is equivalent to that of classifying the texture on 
any bouquet of circles to which it may be retracted^ . For 
each circle in the bouquet, we can choose any element of 
7ri(MP^), so the classification for a bouquet of k circles 
is by a fc-tuple of half-integers. In this case it turns out 
that this is also the answer if we allow the texture at the 
base point to vary. The complication, as we will see in 
the following, is that there can be more than one way to 
tie together the measuring circuits. When 7ri(GSM) is 
non-Abelian this freedom of choice leads to an ambiguity 
in measuring the "charge" of a defect. This is a major 
point of this review and also emphasizes a possibly obvi- 
ous point: the Abelian or non-Abelian nature of defects 
only comes into play when there is more than one defect. 
Making measurements correctly is, as always, the key to 
understanding the physical system. 



C. Three-dimensional Director: RP 

In our two-dimensional example, we insisted that the 
director never point perpendicular to the plane and it 
was then possible to describe the texture in terms of 
idealized configurations where the director was entirely 
planar. This situation is often an accurate description 
of thin samples where the bounding surfaces are treated 
to promote planar alignment. However, in bulk samples, 
the nematic director can point along any of the directions 
in three dimensions with important implications for the 
topology. 

The defects described in this subsection still may be 
captured by a measuring circle in the sample - in a 
three-dimensional nematic, these are line defects. By 



See Mermin's review (Mermin 19791 for a more detailed 



and precise discussion of these properties. 



^ It can be proved that "deformation retractions" of the sample 
onto a smaller space don't change the properties which can be 



probed with the tools of homotopy ( Hatcher 2002 i 



considering a two-dimensional cross-section of a three- 
dimensional nematic, we can discuss measuring circles 
around line defects as circles surrounding points in the 
plane, as in the previous subsection. Thus the sum of the 
measuring circles about two line defects corresponds to 
the process of merging two parallel defect segments into 
one. It is not hard to see that the constructions explained 
in the last subsection have direct analogues for line de- 
fects probed this way, and we leave elaboration on most 
of them to the reader. In later sections, we will probe 
line defects with other measuring surfaces. 

The GSM here is a sphere with antipodal points iden- 
tified (as n and — n are identified in a line field). In this 
instance, the geometry is simple and we are able to see 
in our mind's eye the topology of the GSM. In more gen- 
eral situations, it is useful to have something more sys- 
tematic: indeed, without a systematic approach we can 
never be sure that our intuition is not fooling us. The 
general framework of the Landau theory of phase tran- 
sitions that exploits symmetries and symmetry breaking 
not only provides well-known, systematic tools to our 
problem, but also emphasizes the essential aspect of topo- 
logical defects - a higher symmetry phase surrounded by 
a lower symmetry phase. 

Recall that any rotation in three dimensions is an ele- 
ment of SO {3), and one set of coordinates on this space 
consists of the three Euler angles. Write the rotation ma- 
trices as Ro,^7 = NqM^N^ where the matrices M^ and 



Nj are 


cos 13 sin /3 




cos 7 sin 7 




M^ = 


1 

— sin /? cos 13 


,N,= 


— sin 7 cos 7 
1 


; 






(6) 


and a, 7 


e [0,27r) and /3 


e [0,7r]. 


Away from the 


de- 



feet, the director is a rotation of some fiducial direction, 
say z, so at each point x in space there are three an- 
gles a{x), /3(x), and 7(x), and n = Rq^^z. The uniax- 
ial nematic state has two symmetries. First, because it 
is uniaxial, rotations around n leave the system invari- 
ant (a phase made of cylinders) and so we must identify 
R-Q/37 with Rq^7' so that a rotation around the origi- 
nal z-axis does not change the state. Since we may as 
well set 7' = 0, we first see that we can parameterize all 
directors by only two angles and we only need consider 
the subspace of 50(3) represented by T^^ — Rapo- The 
group of rotations around a single axis is called 5*0(2) 
and this subspace made by Tap is called 50(3)/50(2) 
- elements of 50(3) which are identified with each other 
by 50(2) rotations. Observe that this is a sphere, the 
space of configurations of a unit vector, appropriate for 
a spin in a Heisenberg magnet, and a and /? run over 
the standard azimuthal and polar angles on the sphere, 
respectively. In a nematic, n and — n are further iden- 
tified; they are related by a rotation by n around any 
axis perpendicular to the director n. Again, we can pre- 
rotate the z direction by tt using the diagonal matrix 
P = diag[l,-l,-l], an element of 50(3) with P^ = 1. 



The two element group made of 1 and P is called Z2. 
The nematic symmetry forces us to identify the two ele- 
ments Ta/s and Tq,^P = Ta+Tr.ir-ii- The resulting space 
is SO{3)/H, where H is the group of nematic symme- 
tries (the isotropy subgroup) which includes both 50(2) 
and Z2. Again, the notation SO{3)/H indicates that we 
identify two elements of 50(3), R and R', if they differ 
by any element of the group made oi H ^ {Nq,PNq} 
for all a. Note that PNq,P — N_q, and so the product 
of two elements of H, P™iNci and P'^^N^^ is 

IMqjI' INqs — -^ i^Q2+(-l)'"2Qi- [I) 

While the total number of flips just adds {mi + 7712), 
the rotation angles do not - whether they are added or 
subtracted depends on 7712. Thus the group addition in 
H is not just the group addition of the two different parts 
Z2 and 50(2) separately for to and a. H is called the 
semidirect product of the two groups and, in this case, is 
the point symmetry group Dao, the proper symmetries 
of a right cylinder with constant, circular cross-section. 
In any event, the resulting space SO(3)/H has a name, 
the real projective plane MP , and the identification of 
the spherical angles {9, (f) -^ {tt — 9,(1) + n) allows us 
to view it, just as our intuition suggested, as the sphere 
with diametrically opposed points identified, as shown in 
FIG. [2J However, the way that the angles add ^ may 
not have occurred to the casual observer without this 
analysis'^ . 

This analysis demonstrates that based measuring cir- 
cuits with the topology of a circle in a three-dimensional 
nematic and the one-dimensional defect structures they 

but rather 



may encircle are characterized not by tti 
by 7ri(MP ), which turns out to be a much smaller group. 
Although when we traverse our measuring circuit in the 
sample the final orientation of the director field must be 
the same as the initial one, the final rotation K-afs-y does 
not have to simply be the identity; it can be any sym- 
metry transformation of the nematic, i.e. any element of 
H = Dao- This subgroup consists of two disconnected 
pieces, the rotations Nq and composite rotations PN^. 
Homotopies can alter the precise value of the final rota- 
tion but not which connected component it is in. Conse- 
quently we see that there are two classes of loops in MP^ . 
Returning to the geometric representation of MP in FIG. 
[2] as the sphere with antipodal points identified, we can 
see that these two classes can be identified with loops 
that start and end at the same point on the sphere and 
those that start and end at antipodal points. The pres- 
ence of a defect inside our measuring loop corresponds to 
the latter class of paths. In the previous case, we could 



^ Note that in our discussion of the two dimensional nematic, we 
could have described the GSM as §^ with diameters identified 
and gone through the same construction. The resulting space 
is Rpl, but its topology is identicaJ to that of S^ S^ ^nd RP^ 
are not topologically identical; in particular they have different 
fundamental groups. 




FIG. 2 The GSM for the three-dimensional nematic is the real 
projective plane RP^ . Any measurement around a closed loop 
in the sample can be mapped to a path in RP^ which belongs 
to one of two homotopy classes; trivial (a) and nontrivial (b). 
In other words, 7ri(RP^) = Za. 



simply count the number of brushes to determine that 
7ri(MP^) = iZ, the half-mtegers. In the case of a true 
three-dimensional director we see that the group is actu- 
ally simpler: 7ri(MP^) = l-i since any closed path on MP^ 
can either be smoothly deformed (holding its base point 
fixed) to its base point or to a path that connects an- 
tipodes. We will write %i multiplicatively, so that it has 
the two elements {1,-1}. Thus, no distinction can be 
made any longer between positive and negative strength 
defects; the planar textures Ci/2 and c_i/2 can be con- 
verted into each other by a uniform rotation of all of 
the molecules by tt about some axis in the plane, x say. 
On the GSM this rotation simply corresponds to trans- 
porting the loop n(r) from one side of the equator to 
the other by passing it over the North Pole, as shown 
in FIG. [3] Equally, there are no 'winding 1' defects; the 
combination of circuits around two disclination lines in 
this case is always homotopic to a defect free texture. 
For instance, the combination of the two planar textures 
Ci/2 produces a loop on the GSM that goes all the way 
around the equator and so can be smoothly slid up in 




FIG. 3 The homotopy that takes a -f | disclination line in a 
three-dimensional nematic into a — | disclination line consists 
of a uniform rotation of all the molecules by tt about the x- 
axis. In RP^, this homotopy entails sliding the path from the 
eastern half of the equator (a) over the North Pole (b) to the 
western half of the equator (c) . 



latitude until it shrinks to a point at the North Pole'*. 



^ Strictly this particular homotopy docs not preserve the base 
point, however, it is easy to construct one that does; simply 



The above homotopy from a loop covering the equa- 
tor to a constant loop at the North Pole is actually ob- 
served in cylindrical capillaries treated to give perpen- 
dicular anchoring at the surface. The molecules rotate 
out of the plane so as to align along the cylinder axis in 
the interior - the so called 'escape in the third dimen- 
sion' dCladis and Kleman| |1972[ |Meyer[ |1973| [Wilhams 



and Bouligand, 1974 ). Of course, this escape can proceed 



in two directions: up or down. If, as typically happens, 
the molecules escape up in one section of the capillary 
and down in another then there will be a mismatch in 
between. This gives rise to a different type of defect: 
point defects in three dimensions, known colloquially in 
the field as hedgehogs (Polyakov 1974). 



III. HEDGEHOGS 

Now that we have been introduced to the line defect, 
it behooves us to meet its partner, the hedgehog. Many 
features of the description of line defects carry over to 
the point-like hedgehogs. In this paragraph, we out- 
line what is to come in this section by sketching a few 
"3-dimensional versions" of the pictures described previ- 
ously. First, a measuring loop can link with a defect line, 
but a loop cannot link with a point. Just as we need a 
Gaussian sphere to measure a point charge in electrostat- 
ics, we measure the topological charge of a point defect 
on a measuring sphere (denoted S^) in our sample. Thus 
what we seek to classify are the set of maps without base 
point [S2,IRP^] and the second homotopy group 7r2(MP^). 
With that knowledge in hand, we can consider a piece 
of the sample with multiple point defects by imagining 
a three-dimensional ball with measuring spheres around 
punctures at the hedgehogs. This space retracts to a 
bouquet of spheres, and the "sum" of the textures on 
the measuring spheres in the bouquet is taken to be the 
texture on a sphere which shrink-wraps the entire thing. 
Just as for line defects, this picture relates defect addition 
and a local-to-global construction. 

Let us first study the group structure on 712 (GSM), 
consisting of homotopy classes of based maps from the 
sphere to the GSM. Here again, we fix a base point on 
both the sphere (which corresponds to the base point 
in our sample) and the GSM to define a group addi- 
tion property. To be precise, consider the standard lati- 
tude and longitude coordinates on the globe. We choose 
the North Pole for the standard point on the sphere and 
some convenient point, say Tq, in the GSM. To com- 
bine elements we imagine the following sequence of ho- 
motopies, depicted in FIG.|4J Consider first a single ele- 



hold that point of the loop fixed on the equator, slide the rest 
over the North Pole and shrink it into the base point. In this 
case too, [S-'jRP^] has the same elements as 7ri(IRP^). In general 
[S^jGSM] is equal to the underlying set of the group 7ri(GSM) 
whenever the latter is Abelian. 



ment gi G 7r2(GSM). We can smoothly deform the map 
by combing out a neighborhood of the North Pole into 
a polar cap which maps entirely to Tq. We can keep 
smoothly deforming the whole map into a smaller and 
smaller patch of the sphere, until we have a small island 
surrounded by an ocean of points all mapping to Tq. 
The texture arising on the island from a standard -|-1 
hedgehog is a skyrmion (Miihlbauer et al.\ 2009 Rofiler, 



Bogdanov, and Pfleiderer 2006 ) and the director rotates 



from To at the coast all the way to the antipodal di- 
rection Tq in the interior, covering every orientation in 
between. This new map is homotopic to the original map 
which corresponded to gi and so, from the point of view 
of 7r2(GSM), the new map is the same element gi - re- 
member, the rate at which we traverse the GSM is not im- 
portant, only the places visited. We can make the same 
smooth maneuver on a second element 92 G 7r2(GSM), 
making a second island in the same ocean. To combine 
the two elements, we simply put island one (gi) and is- 
land two (32) in the same ocean without overlapping the 
islands. This works because the boundaries of these is- 
lands are all mapped to Tq. We can now deform this new 
map smoothly as we see fit, filling the ocean back up with 
land. In doing so we have combined the two elements. 
Since we can move the islands around before filling in the 
ocean, there is no way to order the group addition and 
so 7r2(GSM) and it follows that tt2 is necessarily Abelian. 
As in the case of tti, the identity element is again the 
equivalence class of the uniform texture that maps the 
whole sphere to Tg. 

Finally, the inverse of any element can be found by 
deforming the original map into an island and then flip- 
ping the island over on the sphere. Why docs this give 
an inverse? We can connect each point on the Western 
hemisphere to its reflection through the plane including 
the prime meridian in the Eastern hemisphere. We can 
choose the value of the map on each of these lines in- 
side the two sphere to take the value at the endpoints 
- identical by construction. It follows that we can fill 
the region inside the sphere with a smooth texture and 
there can be no defects inside. Thus the island and its 
mirror add to zero, shown in FIG. [5] Alternatively, we 
can take the map f{9,(t>) i-s> GSM and create the inverse 

Let's relate the group structure of 7r2(GSM) to the ad- 
dition of defects in terms of measuring surfaces: imagine 
(as in FIG.l4| that we measure two different regions, each 
with a sphere. As for loops, in order to add the sphere 
measurements, they must be in a bouquet configuration. 
We can achieve this by attaching a tubelike tether to each 
sphere ending at our base point. Beware! The charges 
assigned to hedgehogs (actually, assigned to measuring 
spheres) will depend on the choice of tethers, particu- 
larly if there are disclinations in the system. The sum 
of the textures on these spheres is equivalent to the tex- 
ture on a large sphere which has been deformed so that 
it snugly wraps around all of the original spheres. We 
encourage readers to relate this to the earlier picture of 




(e) r'^ 




FIG. 4 How do you combine two hedgehog charges (a)? The key hes in based hoinotopy theory. First, we comb the texture 
on the measuring surface so that the defect is confined to a small island within an ocean of uniform director field, which we 
identify as the base point (b). In order to compare the relative charges of two hedgehogs, their respective measuring surfaces 
must originate from the same base point, forming the bouquet of spheres described in the text (c). A new measuring surface, 
a sphere which has been shrink-wrapped around the bouquet of spheres, now contains both islands (d) and may be deformed 
into a single sphere containing the sum of both hedgehogs (e). 



islands on a single sphere - each sphere in the bouquet 
becomes an island, and the base point of the bouquet is 
blown up to form the ocean. If no disclination lines are 
present, i.e. ■7ri(GSM) — {1}, then the tethers can al- 
ways be unentangled since you "can't lasso a basketball" 
( |Colemannl988) . 

Though we have a definition and interpretation of 
7r2(MP^), we haven't computed it yet. In order to do 
so, we generalize the idea of counting brushes in the two- 
dimensional nematic. Each brush is a place where the 
director points in one of two directions - parallel to the 
polarizer or analyzer - and the charge is one quarter the 
number of brushes. Because of the cross-polarizers, we 
are forced to count both the polarizer and analyzer direc- 
tions. However, were we given an explicit map we could, 
instead, just count the number of times the map from 
S^ went to any single element of MP^ (with sign, as in 



Eq. (151)), an integer known as the degree of the map. 
Then the charge of the defect would be half the degree. 
We could try to do the same thing for maps from S^ to 
MP^, but the ambiguity of n and — n makes this prob- 
lematic^. Each point of RP can be identified with two 
antipodal points on §^ : we can start at the base point To 
and choose one of the two equivalent points on §^ . Given 
a choice, we can do this for all points in a small neigh- 
borhood consistently so that we locally generate a new 
smooth map to S^. The only hitch might be that when 
taking a long path around MP we may find that we are 



^ The reader may worry that the same problem plagues us for 
7ri(]Rpl). It does not: MP^ is oricntable where MP^ is not, which 
is the real issue here. 
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wraps around the sphere, also known as the degree or 
hedgehog charge 



FIG. 5 The inverse of an element consists of taking the mirror 
image of its island and placing this new island on the sphere. 
The sum of these two elements must add so that the sphere 
contains no defects. In the interior of the sphere, the texture 
is constant along chords connecting identical points on the 
original island and its mirror. Because no defects are intro- 
duced, the entire texture on the sphere may now be smoothly 
deformed to the base point. 



hung up on a non-contractible path in the GSM when we 
return to the same point^. However, this is impossible 
with maps from §^. Fortunately, maps from §^ to any 
space X have a special property owing to the topology 
of the sphere. Consider a closed loop F on S^ and the 
loop it produces on X and watch how the latter changes 
as we smoothly deform F. Because F can be shrunk to 
a point on S^, its image on X can also be shrunk to a 
point and to one value in X. For the case of KP this is 
very convenient: the image of any loop F on §^ must be 
homotopic to the identity map on MP , so we never get 
hung up. Once again, we cannot lasso a basketball and 
so we can always "lift" the map from MP to S^ globally 
and turn it into a map from §^ to S^. Therefore, after 
fixing the lift of Tq (for definitiveness, we'll take it to 
lie in the Northern hemisphere), 7r2(MP ) is the same as 
7r2(§^), which we can calculate. 

Each map gives us a unit vector n{9, (p) at each point of 
the measuring sphere. Like the winding number, there is 
an integral that measures the number of times our map 



This is the distinction between the identity and — 1 in 7ri( 



')■■ 



1 

47r 



§2 



n • [den x cJ^n] 



(8) 



Remarkably, it is an integer: the integrand is just the Ja- 
cobian of the map from S^ to S^ and so it counts the area 
swept out on the target sphere. Dividing by 47r simply 
gives us the number of times we visit each point on the 
target. Moreover, because we do not take the absolute 
value of the Jacobian, we measure positive and negative 
area and thus will properly count the analog of increasing 
and decreasing as needed in ([5| fo r the winding number. 
See, for instance, (Kamien 2002) for the derivation of 



this Jacobian and its connection to Gaussian curvature. 
Therefore, we see that maps from §^ to S^ are classified 
by a degree and that this degree can be any integer, pos- 
itive or negative. Thus 7r2(MP ) = 7r2(S^) = Z and there 
are an infinite number of topologically distinct point de- 
fects in nematic liquid crystals. Representatives of each 
homotopy class are given by the maps 

nrf(6', (f>) = [ain{9) cos(d(/)), sin(6') sin(d0), cos(6')] , (9) 

which exhibit a d-fold winding on the equator and may 
be thought of like the textures c™ of equation (l2| as 
providing an idealized description of the director field 
near a hedgehog. 

There remains just one point to straighten out: the 
question of the choice of lift. This will also tell us about 
the difference between the free homotopy classes [S^ , MP ] 
and the second homotopy group 7r2(MP^). Suppose we 
had instead chosen to lift the base point Tq into the 
Southern hemisphere of S^. What would be different? 
We would still be able to identify each element of 7r2 (MP ) 
with one in 7r2(S^) and define a degree d through the for- 
mula Q. However, when we would have used the vector 
n{6, (f)) to specify our texture we would now use — n(0, (f)). 
Since equation (Is]) is odd in n the degree we now mea- 
sure would be —d instead of d. This is not to say that the 
textures labelled by d and —d are the same; the identifica- 
tion we made between 7r2(MP ) and 7r2(S^) is one-to-one, 
it's just that there are two identifications that we can 
make, differing in whether we lift the base point to the 
Northern or Southern hemisphere, and these are not the 
same^. So long as we make our choice of lift consistently, 
hedgehogs with charge d and —d are distinct, but it is up 
to us which sign we associate to which hedgehog®. 



when paths in RP^ are lifted to paths in S^, paths in the equiv- 
alence class of the identity are closed, and paths in the class —1 
are not. 



The reader might be concerned that the same ambiguity arises 
for 7ri(Rpi). However the two choices of lift from KP^ to §1 differ 
by a TT rotation, i.e., to replacing ft(<t') with ft{<l>) + ''r in equation 
(J5l, which does not change the winding number w. In this case 
both lifts yield the same identification of 7ri(RP ) with 7ri(S"'^). 
The prevailing convention in the literature is to consider the 
purely radial texture, equation ^ with d = 1, to have posi- 
tive charge, the choice of lift being that the director field points 
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However, under free homotopy maps from the sphere 
to MP labeled by d and —d do become equivalent - take 
the texture of (|9| and rotate uniformly by n around the 
y-axis, then take n to — n to change d to —d. It fol- 
lows that [S^jMP'^] is the set of nonnegative integers, 
which is not the set of elements in tt2 (MP ) ! Thus even 
with the same spherical measuring surface, the classifi- 
cation of hedgehog defects can be different depending on 
how we make the lift from MP^ and S^. As long as we 
don't try to extend the local measurements to a global 
texture, the "local-to-global problem", we need not be 
careful. However, if we specify free homotopy classes on 
spheres about defects and then attempt to fill the remain- 
ing space, we run into ambiguities - trouble! Were we to 
probe each hedgehog individually by calculating the de- 
gree on a small sphere around it, we could only measure 
its charge up to a sign and we are assigning it to one of 
the maps [§^ , MP ] . The ease with which we do this ex- 
tracts a penalty - we can no longer add defects since we 
have lost the group structure of 7r2(MP^) without the use 
of a base point. There is no free lunch, regrettably. For 
instance, a pair of spheres with \d\ = 1 on both might 
induce either \d\ =0 or 2 - we can't know the relative 
signs unless we have a common reference coming from a 
fixed base point. 



To understand the global structure it is especially il- 
lustrative to view the sample as a bouquet of spheres 
attached at a single point, the base point of the sam- 
ple, that maps to the base point of the GSM. If we have 
k point defects then the texture becomes a map from 
a bouquet of k spheres to MP . Were we to measure 
the "charge" on each sphere using the base points, we 
would end up with a fc-tuple of elements of 7r2(MP^) 
labelled by their degree, {di,d2, ■ ■ ■ ,dk). There is al- 
ways a global ambiguity when we lift to S^ since the 
GSM base point can be taken to lie in either the north- 
ern or southern hemisphere. If we now choose to con- 
sider the free homotopy classes of samples with k defects, 
we identify {di,d2, ■ ■ ■ , dk) with (— di, —d2, ■ ■ ■ , —dk), not 
(±(ii, ±(i2, • ■ • 7 irffc) In other words, it is possible to go 
from a based measurement of a single defect to a based 
measurement of a collection of defects to an unbased mea- 
surement of many defects. It is not, however, possible to 
make an unbased measurement of many defects from un- 
based measurements of single defects. This is the crux of 
the local-to-global problem and is why local information 
is not always enough to understand the topology of the 
whole texture. 



IV. DISCLINATION LOOPS 

We have now described disclination lines and hedge- 
hogs separately, but in any typical texture both will be 
present simultaneously and we are forced to think about 
how they interact with each other. When we discussed 
defect lines previously, we merely generalized point de- 
fects into a third direction and classified them by mea- 
suring circuits which ensnare the defect lines. However, 
these measurements could not determine whether or not 
the disclination itself is a loop or a knot, or linked with 
other disclinations. Indeed, a great part of the appeal of 
colloidal systems is precisely that they provide a natural 
setting for illustrating the topological interplay between 
disclinations and hedgehogs and enable the braiding and 
tangling of the disclination lines (Tkalec, et al. 2011). If 
the colloid is treated to promote radial anchoring of the 
molecules at its surface then it will appear like a unit 
strength hedgehog and this topological charge will have 
to be compensated by a companion defect in the liquid 
crystal so as to satisfy the global boundary conditions. 
Another way of saying this is that the homotopy class of 
the texture on a measuring sphere which surrounds just 
the colloid differs from the one on a much larger measur- 
ing sphere on which, say, the texture is determined by the 
global boundary conditions. This implies the existence 
of another defect set. One way of achieving this balance 
is to place an opposite strength point defect close to the 



colloid so that together they form a dipole pair (Poulin, 



et al. 1997) . Multiple dipolar colloids interact with each 
other elastically through the distortions they produce in 
the liquid crystal and can as a consequence a ssemble into 
cha ins or two-dimensional colloi dal crystals (Musevic, et 



al. 



2006 



Skarabot, et al. 2007). 



A separate possibility is to have the charge of the 
colloid compensated by a disclination loop, encircling 



the particle in a 'Saturn ring' configuration (Terentjev 
1995). These too interac t elasticallv and aggre gate to 



form chains and crystals (Skarabot, et al. 20081. How- 



ever, they also allow for a variety of intriguing entan- 
gled structures where a single disclination loop wraps 
around two or more colloids, balancing their collective 
charge (Guzmarij^^e^ al.\ 2003 Ravnik, et al. 2007 Ravnik 



and Zumer, 2009). This naturally begs the question of 



how to describe the topological properties of such discli- 
nation loops, which may be probed both as line and point 
defects. Although this makes sense from the perspective 
of conserving total charge, hedgehog charge, it turns out, 
is not a homotopy invariant of a the full measuring sur- 
face of a disclination loop - a torus. 



A. Measuring with Spheres 



outwards. The base point on the measuring sphere would be the 
North Pole and the base point in RP is the vertical direction, 
which is then chosen to lift to the North Pole in §^. 



Consider a based measuring sphere in our sample 
which doesn't surround any defect. The nematic texture 
can be deformed so that the sphere has one island with a 
-f 1 charge and a mirror island with charge -1. As shown 
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FIG. 6 A sphere with a plus and minus island has zero net hedgehog charge (a). The two islands can be separated to create 
a plus {S+) and minus (S_) pair (b). As the minus hedgehog is carried around a disclination line its orientation reverses as 
the texture on the tether between the two dumbbells connects antipodal points of RP^ and the natural final measuring circuit 
{S'_) records the opposite hedgehog charge (c). If the two hedgehogs are now recombined and the tether is ignored their total 
charge will be ±2! (d) 
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in FIG. [6J we can imagine distorting the sphere into a 
dunibbeU shape with one island on each end. As we pinch 
off the tether connecting the two ends, each island comes 
to reside on its own sphere, 5*+ and S-, surrounding a 
plus and minus hedgehog, respectively, which we real- 
ize were necessarily created during the process. Suppose 
that elsewhere in the sample there is a 1/2 disclination 
line. What happens if we drag the -1 hedgehog around 
that disclination? More precisely, we keep the texture 
fixed on S+ so that the degree integral is unchanged but 
begin deforming S- so that it wraps once around the 
disclination line and returns to its original position, leav- 
ing a tail tethered around the disclination line. Now, the 
most natural sphere around the -1 hedgehog is one that 
doesn't have the tether, S'_, and is the one which we 
would likely use to measure the charge locally. What's 
the relation between the textures on S- and S'_7 It can 
be difficult to visualize this operation, which involves a 
homotopy of a three-dimensional line field, but if we re- 
turn to the original undeformed texture, the texture on 
S'_ is in fact equivalent to the one on S- except that 
we then have attached a loop to it which goes around 
the disclination line the "opposite direction" to the base 
point. Roughly speaking S- and S'_ clasp the hedgehog 
from different sides of the line defect. Let us call this set 
of operations maneuver X. Such a tether picks up the tt 
rotation of the disclination and we are forced to change 
n to — n on most of the sphere, in particular over the 
island, when we lift the texture on S'_ from MP to S^. 
Thus when we calculate the degree on S'_ via integration, 
d will go to —d\ We have reversed the degree of the map 
on the sphere around the hedgehog by "moving it around 
a disclination." Note that we had to carefully compare 
two distinct measuring spheres around the hedgehog in 
question in order to make sense of this. We started with 
total charge around the two hedgehogs and ended with 
total charge ±2 because we changed the measuring sur- 
faces - just as the total charge in a Gaussian measuring 
box cannot change, neither can the charge in fixed mea- 
suring circuit. 

We may also interpret maneuver X from the point of 
view of "islands on the globe." Imagine a disclination 
loop that approaches the globe. We may poke an island 
on the sphere through this loop, and while performing 
this process the loop leaves an image in the shape of an 
atoll in the Tg-ocean. Though the value in MP^ can be 
the same inside and outside the atoll, when we lift to S^ 
we must change the sign of n when we cross the atoll, 
as this is equivalent to wrapping around the disclination 
loop. If the atoll is only surrounding a region of open 
ocean, this is not a problem - we can shrink it to a point 
and effectively make the whole region of "wrong" ocean 
of Tq disappear. But given an island carrying a non- 
trivial hedgehog charge, when we lift to S^ we will be 
forced to match the island's coastline to a lagoon of Tq. 
No problem, just lift the island to — n instead of n. Is 
there a singularity of the texture on the atoll? No! For 
concreteness take Tq to be the North Pole and so Tq is 



the South Pole. Topologically, the atoll is an annulus. 
The inner ring of the annulus points South and the outer 
ring points North. It is no problem to have the direc- 
tor smoothly rotate along the radius of the annulus from 
South to North. This texture will not contribute to the 
degree. By stretching the measuring surface this becomes 
precisely the picture in FIG. |6] where the atoll becomes 
the tube-like tether. Again, note that this action of tti 
on 7r2 preserved the island's class in the set of unbased 
maps [S2,MP^], but not in 7r2(MP^). 



B. Measuring with Tori 




fA^^ 



FIG. 7 The measuring surface T^ has two cycles, either of 
which may contain nontrivial winding. For the unlinked 
disclination loop shown the winding is non-trivial around the 
meridian and trivial around the equator. 



Hedgehog charge is a homotopy invariant of a sphere, 
the natural measuring surface for a point defect. But 
the natural measuring surface for a disclination loop is 
not a sphere but a torus, T^, a thin tube sheathing the 
singular line. The classification of disclination loops can 
therefore be based on the homotopy invariants of the tex- 
ture on this torus, that is on maps from T^ to MP^. There 
are multiple measures of these maps, two corresponding 
to elements of 7ri(IRP^) measured on the cycles of the 
torus and a more refined quantity, a global defect in- 



dex (Bechluft-Sachs and Hien 11999 Janich 1987 Nakan- 



ishi, Hayashi, and Moril 1988[ ) which captures some as- 
pect of hedgehog charge. A torus surrounding a simple, 
circular disclination loop, shown in FIG. [Tja, has two 
cycles: a meridional loop that goes around the disclina- 
tion line which always measures the non-trivial element 
of 7ri(MP^) and the loop that follows the contour of the 
disclination line along the longitude of the torus records 
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another element of 7ri(IRP^)^. If this is also non-trivial 
then it means that our disclination loop itself goes around 
another defect - linked loops! - while if it is the trivial 
element then our disclination loop is either isolated, or 
linked an even number of times. For now let us assume 
that this element is trivial and that our disclination loop 
is unlinked. 

What is the homotopy classification of this subset of 
textures on a torus? One naive guess would be that we 
would have an additional choice of integer n, which we 
might get by repeatedly merging hedgehogs into the line 
defect, in a way analogous to our way of merging two 
hedgehogs together (via a bouquet construction). How- 
ever, we can show that adding two hedgehogs to our 
disclination is the same as not adding any by using ma- 
neuver X, as in FIG.|8] Let us have a texture on the torus 
called /; we will now show that /' arising from adding 
two -t-1 islands to the torus by merging in a sphere that 
carries +2 hedgehog charge is homotopic to /. Move one 
of the -|-1 islands on /' on a loop around the meridian. 
As we pull the island around a circle of longitude it accu- 
mulates a concentric series of borders around the island 
which rotate the coastline by tt. As before this is the 
action of -1 e 7ri(MP^) on +1 e 7r2(MP^). This island 
is now homotopic to — 1 e 7r2(MP^), as we argued before. 
Now we have a ±1 pair of islands on the torus which 
we may cancel against each other just as we could for a 
pair of such islands on a sphere. This results in a texture 
on the torus which is homotopic to the one we started 
with. From another point of view, performing maneu- 
ver X does not require any discontinuous changes on the 
torus - the hedgehog moves through the handle with- 
out piercing the measuring surface. As a result, we can 
smoothly transform the texture on T^ while changing the 
hedgehog charge of the surrounding space by ±2. Hence, 
we can classify tori surrounding disclination loops as ei- 
ther carrying even or odd hedgehog charge. The reader 
might ask how we can be sure that there are no moves 
that can change this charge by ±1, an issue that we will 
explain in the next section. 

In what follows we will state some interesting results 
that fill out most of the rest of the story about discli- 
nation loops. The complete classification of textures on 
tori labels textures with one of the four pairs (a, b) where 
a,b = ±1 are the homotopy classes on the meridian and 
longitude cycles, respectively. For the three cases where 
a and b are not both +1, it turns out there are two 
subclasses of textures within those that are labeled with 
(a, 5), which correspond to even and odd numbers of -1-1 
islands on the torus, just as for the case (—1, +1) we ex- 
plained above. When a and b are both -1-1, there are an 
infinite number of classes - without a disclination loop 






® The reader may worry about which longitudinal path to follow if 
the disclination loop has a more complicated shape or is knotted. 
A canonical choice is provided by a ny path that has zero linking 
number with the defect loop itself I Janich 1987 1. 



EAh 



FIG. 8 Two -1-1 islands on a torus surrounding a disclination 
loop, representing the addition of two hedgehogs (a). Moving 
one of the islands around the meridian of the torus reverses its 
orientation, turning it into a charge —1 island (b). If the two 
islands are then merged they carry zero net charge, illustrat- 
ing the homotopy between disclination loops with hedgehog 
charge n and n ± 2 (c) . 
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in the torus, we cannot cancel out pairs of islands any- 
more^'^. We therefore can add a subscript either in {e, o} 
or Z to this ordered pair to complete the classification. 

This peculiar set of homotopy classes has some in- 
teresting additional structure: we may do much more 
with tori than just merge them with spheres, as we did 
above. We present below simply one example, which cor- 
responds to merging two disclination loops (i.e. tori with 
a = — 1 ) "sid e-to-side" (Bechluft-Sachs and Hien 1999 
Janich 1987[ ). Unfortunately, a precise explication of 
the other group which emerges when we merge two un- 
linked tori (i.e. with 5=1) "top-to-bottom" (Nakanishi, 
Hayashi, and Mori 1988[ ) is just barely outside the scope 
of this paper. Given two tori with textures such that 
a = —1 on both, we cut them along a meridional circle, 
resulting in two cylinders and then reglue them so that 
we have a single torus. This results in a Z4 group struc- 
ture on the homotopy classes of such textures on a torus, 
where the mapping is [0] = (— l,l)e, [1] = (— 1,— l)o, 
[2] — (—1,1)0, and [3] — (— 1,— l)e, where we write Z4 
additively so that [?7i]-|-[n] = [m + n mod 4]. It's natural 
for (—1, l)e to be the identity element; it is both unlinked 
and may contain no hedgehogs, which makes adding it a 
bit like adding a "constant" segment of disclination line. 

We summarize these and the results of the last three 
sections in Table IIV.BI 



V. BIAXIAL NEMATICS AND THE ODD HEDGEHOG 

In the last section, we argued that by a smooth de- 
formation, a map from T^ to EP^ could absorb hedgehog 
charge in pairs so that the hedgehog charge could only be 
even or odd. Here, we probe this further and use the in- 
sight provided by decorating the uniaxial textures with a 
small amount of biaxial order. This biaxial point of view 
both highlights the underlying topology of the uniaxial 
phase and clarifies the way in which the even and odd 
classes of uniaxial disclination loops are distinct. 

Recall in our discussion of the nematic we discovered, 
via some matrix algebra, that the ground state manifold 
was MP . This was appropriate for a uniaxial nematic, 
which had rotational symmetry around its long axis ~ the 
rotational symmetry responsible for the 5*0(2) factor in 
H. A biaxial nematic has a lower symmetry, that of a 
brick or rectangular cuboid with three unequal lengths. 
Not just a mathematical construct, biaxial liquid crys- 
talline phases have been known for many decades. In the 
past few years, discovery of thermotropic biaxial phases 
has renewed interest in their defect structures. Moreover, 



^^ Free and based homotopy classes have the same classification 
for the first three classes where a, b are not both trivial, since 
multiplying by —1 doesn't change evenness or oddness of the 
hedgehog charge on the islands. For the (-1-1, -1-1) case, the based 
classification is all integers, and the free classification is only 
nonnegative integers, just as for textures on a sphere. 



chirality and biaxiality are intimately connected (Har- 
ris, Kamien, and Lubensky} |1999| [Priest and Lubensky[ 
1974[), and stu dies of blue phases (Dupuis, Marenduzzo, 
and Yeomans[ [2i305[ [Grebel, Hornreich, and Shtrikman 



1983 Wright and Mermin 19891 often utilize a biaxial 



description. 

Returning to the notation and discussion in Section 
|II.C[ we start with the biaxial molecule with long axis 
along z and a second axis along x (the third axis is along 
y and all three are distinguishable - "triaxial nematic" 
might be more apt). The symmetry of the brick involves 
only three discrete rotations of tt around each of x, y, 
and z. Ignoring those symmetries for the moment, we 
see that the original rotation matrices Rq^-^ represent 
an arbitrary rotation of the brick. With the symmetries, 
we must now identify 7 with 7 ± tt, so that the isotropy 
subgroup Hh becomes 



iffc = {l,P,N,,PNJ, 



(10) 



where the addition of elements (multiplication of the ma- 
trices) is as before. Note that PN^r = M^r, so that a ro- 
tation around z of tt followed by a similar rotation about 
X yields a tt rotation around y. For all three of these ele- 
ments, we are never sure whether these rotations are by tt 
or — TT. Of course, the reader might think that these lead 
to the same group elements, which they do. However, 
a loop in 50(3) that starts at rotation and ends at a 
27r rotation is not contractible and leads to a non-trivial 
element of 7ri[50(3)/iJb]. This fact, one might recall, is 
often demonstrated in a class on quantum mechanics by 
someone who takes off their belt or holds a filled coffee 
cup with one hand and performs an elegant gyration of 
their arm^^. It is why spinors must change sign under 
rotations by 27r and why the spin and statistics of par- 
ticles are interrelated. Now that we have introduced the 
notion of a lift when discussing the hedgehog charge, it 
is simplest to demonstrate this fact with yet another lift. 
The Pauli matrices: 



o-i = 



1 

1 



(T2 = 



-i 
i 



0-3 



1 
-1 



, (11) 



Sijl as follows from their com- 



satisfy (t,(Tj = ieijkcrk 
mutators and anti-commutators. This can be used to 
form a simple way to parameterize rotations in three- 
dimensions. Write any vector x — [xi,X2,X3] as the 
matrix x = x^tTk where we employ the summation 



convention over repeated indices. Then x 



1 is 



the unit matrix times the squared length of the vec- 
tor. Moreover, if U is any 2x2 unitary matrix, then 
if we define the similar matrix x' = U^^xU, we find that 
(x') = U^ja^plU — x^ so under a unitary transforma- 
tion the magnitude of the matrix is unchanged, precisely 



Some readers may have even been subjected to "spinor span- 
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Physical system 



Measuring 
GSM Circuit 



Classification / Notes 



Point defects in Schlieren texture 
Disclination lines in 3D nematics 
Point defects in 3D nematics 

Disclination loops in 3D nematics 



§1 

§1 X §1 



|Z, when based, can be added 



II.A 



II.C 



Z2, when based, can be multiplied 

N when free, Z and may be added when based 

attaching a tether which runs over — 1 G 7ri(RP^) takes d to —d 

(a, b)p, where a,b = ±1 and p = even, odd if a , b are not both 1, if 
a = b 



For based textures, 
IV 



1, then p G Z if based, or N if free 



IV 



TABLE I A summary of classifications of textures on various measuring surfaces, their physical relevance, and how 
they may be combined. If free or based is not specified, the underlying sets for the classifications are the same. 



what we need for rotations. We finally note that the ma- 
trix XJ{0) = 1 — i-^OkCk generates infinitesimal rotations 
of X since to linear order in 9 (note that these angles are 
not the same as the Euler angles a, /3, and 7): 



U^xU 



i^OkO-k 



i-9-i(Ti 



iXit/kcrk(^i - iXittkO-iCTk 



[Xj+X^eke^k^(T^+0{e'^), 



0{9^) 



(12) 



resuhing in precisely the expression for the rotated vec- 
tor. The fuU rotation is^^ 



U _ g-J0fcCTfc/2 



COS 



. fm\Ok 

3m — -^CTk. 

\^J\e\ 



(13) 



Thus, each rotation in 5*0(3) can be identified with a 
2x2 unitary matrix with unit determinant, the group 
SU{2). However, both U and — U generate the same 
rotation or, in other words, 9 and —9 generate the same 
rotation, and thus, like the lift of KP to S^, we have a 
sign ambiguity. Starting at 9 = 0, we can move along 9i, 
holding 02 = ^3 = 0. In this case as 9i goes from to 27r, 
U(0,0,0) = -U(27r,0,0). We end up at the identical 
place in 5*0(3), but at a different place in SU{2). It 
follows that this loop is not contractible to a point and so 
the map from the measuring circuit to 5*0(3) which winds 
by 27r is one of the non-trivial elements of 7ri[50(3)], 
traditionally denoted as —1. 



^^ To establish this exponential formula, consider the unitary ma- 
trix S which diagonalizes the Hermitian matrix t = Sj.(Tfc/2 = 
SAS^ where A is the diagonal matrix of eigenvalues. Then T = 
exp{it} = exp{iSASt} = Sexp{jA}St. Since Tr A = Tr t = 0, 
A = diag[A, —A], and it follows that 



e-'^ 



St=cos{A)l + i''^SASt, 
A 



cos(A)l + 2sin(A) — . 
A 



Dividing out by Hi, yields additional defects akin to 
the stable nematic defects of charge 1/2. In this case, 
however, there are tt rotations around each of i, y, and 
z. Lifting those to 5C/(2), we find that they correspond 
to ±i<Ti, ±i<T2 and ±10-3, respectively, where the sign 
ambiguity corresponds to the same lifting ambiguity as 
above. These group elements have a nasty habit - they 
don't commute, so a rotation of tt around x followed by 
a rotation of tt around y, followed by a rotation of — tt 
around x gives [«ci]~ [1(^2] i*""!] = 'icriC2Ci = —io'2- 
This means that the combination of paths on 50(3) that 
begin and end at one of the elements of H\j do not com- 
mute and TTi[SO{3)/Hh\ is non-Abelian. 

What about point defects in a biaxial nematic? As 
it turns out, there are no point defects in the biaxial 
system, so n2[SO{3)/Hb] = {1} the trivial group of one 
element. To see this only requires a little more work 
and a theorem that we will not prove but is covered in 



standard texts ( Hatcher 2002 ), including the gold stan- 
dard by Mermin (Mermin 1979). First we ask what does 



SU{2) look like as a manifold? Fortuitously, we have the 
answer in the above algebra. Note that we can consider 
a general matrix depending on the four real parameters 

i= {xo,x) eR'^: 



V = a;ol + ixkcrk 



xq + ixs ixi + X2 
ixi — X2 Xq — ix-j, 



(14) 



so detV — 



is the squared length of the vector 



X in four dimensions. Since all elements U S SU{2) are 
of this form and detU — 1, we see that SU{2) can be 
identified with the three-dimensional sphere S"^. So what 
do maps from S^ to S'^ look like? They are all the same 

- just as 7ri(S^) vanishes because we can shrink every 
loop on the sphere to a point, so too can we shrink every 
two-sphere on a three-sphere to a point. Thus 7r2(§'^) = 
TT2[SU{2)] ^ {1}. To calculate 7r2[5*0[3]/7?f,], we need 
only use the theorem that 7r2(G) = 7r2 (G/iJ) if i/ is a 
discrete group^'^. Note that we saw this with S^ and MP 

- they both shared the same ■n2- 



Since — A = det A = det t 



LJ— , the identity follows. 



For the cognoscenti recall that there is a long exact sequence of 
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FIG. 9 Schematic of biaxial order on a colloid with radial 
anchoring conditions for the long axis. The two shorter axes 
undergo a 27r rotation at the North and South Poles and the 
particle is threaded by a —1 disclination so that a Saturn ring 
defect accompanying the colloid is linked with this biaxial 
disclination. 



Precisely because the biaxial nematic does not support 
point defects, it serves to illuminate our previous discus- 
sion of disclination loops and their relation to hedgehogs. 
To see what we have gained it suffices to consider the sim- 
plest example: a disclination loop surrounding a single 
colloid in a Saturn ring configuration, illustrated in FIG. 
[9| As usual the meridional cycle of a torus sheathing the 
disclination records the non-trivial element of 7ri(MP^) 
and here, since the disclination is not linked with an- 
other, the longitudinal cycle records the trivial element. 
However, this cycle is still interesting. Note, in particu- 
lar, that the director field on the inside of the torus (the 
part closest to the colloid) undergoes a 27r rotation as we 
traverse this longitudinal cycle. Thus, if we were to deco- 
rate the uniaxial texture with a perpendicular short axis 
to produce a biaxial nematic we would find that this cy- 
cle corresponded to a non-contractible loop in S0{3)/Hb, 
recording the element —1 of 7ri(50(3)/ifb). The Saturn 
ring was linked with another disclination after all! An all 
but invisible one only revealed by adding a small amount 
of biaxiality. This additional biaxial defect is not present 



homotopy groups 

...-^7T2{H)- 



■7r2{G)^TV2{G/H)^7ri{H) 



for all defect loops: it is associated, as anticipated at the 
end of the last section, with a disclination loop that car- 
ries odd hedgehog charge, of which the Saturn ring is the 
simplest example. 

It is useful to explore this connection a little further. 
Suppose we remove the colloid, leaving behind just the 
disclination loop and the texture on the torus surround- 
ing it. Could we fill things back in differently, say without 
the hedgehog? Note that the configuration on the inside 
of the torus is just the same as on a cylindrical capillary 
with perpendicular anchoring. Thus to avoid a singular- 
ity the director will have to escape in the third dimension, 
just as it does in a capillary. However, if we have biaxial 
order then the singularity cannot be avoided; although 
the long axis is well defined along the capillary axis the 
two shorter axes are not and we have a —1 disclination^*. 
Let us assume that the texture is uniform at large dis- 
tances so that this defect too closes up to form a loop. 
Technically, we are considering textures with defects on 
S'^, the usual three-dimensional Euclidean space with a 
point at infinity added. The biaxial —1 disclination can 
also be surrounded by a torus, whose meridian records a 
2tt winding and longitude records the tt winding of our 
original 1/2 loop with which it is linked. Now by escaping 
in the third dimension the long axis can be made regular 
on any local section of the —1 defect, but this cannot 
be extended globally since the long axis is required to 
undergo a tt rotation along any path encircling the 1/2 
disclination. This tt rotation reverses the orientation of 
the director and converts an initial escape 'up' into an 
escape 'down' so that there must be a mismatch some- 
where. Just as in the cylindrical capillary this mismatch 
marks the location of a point defect in the uniaxial ne- 
matic - the odd hedgehog. 



VI. CONCLUSION 

In the hope of avoiding to provide only a poor imita- 
tion of the many excellent reviews of the homotopy the- 
ory of defects already available, e.g. (Kurik and Lavren- 
tovich 



il988, Mermin, 1979, Michel 



1980: Trebin 19821, 



we have focused our discussion on a single physical sys- 
tem, nematic liquid crystal colloids, rather than provide 
a general survey and have eschewed, as far as possible, 
all of the formalities of algebraic topology. 

A recurring issue is the necessity of a base point to in- 
duce the group structure on the set of homotopy classes. 
From a physical perspective the base point is really a 
fiction, a convenience introduced for its useful compu- 
tational attributes rather than because it conveys any 
deep significance. As mentioned in Section |TT] the choice 
of base point for this purpose is entirely arbitrary so that 



and if H is discrete then ■K2{H) = tti{H) 
follows that 7r2(G) = ■n2{G/H). 



{1}, from which it 



^* This transferral of winding from one a xis to a perpendicular o ne 
is the celebrated Mermin-Ho relation (Mermin and Ho 1976 i. 
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any choice we make should yield the same results as any 
other. Fortunately, it is a simple exercise, proved in all 
of the standard reviews (Mermin 19791, to show that 
the homotopy groups 7ri:(GSM, Tq) and 7rfc(GSM, Tq), 
defined with base points Tg and Tq, respectively, are iso- 
morphic. The physical properties of defects should be in- 
dependent of any choice of base point that we may make 
and the isomorphism between homotopy groups with dif- 
ferent base points assures this. More importantly, how- 
ever, a continuous deformation of a liquid crystal texture 
need not, in general, hold any point on any measuring 
circuit fixed while all others can freely vary. No physical 
distinction should be ascribed to defects that are non- 
homotopic when a base point is used, but become homo- 
topic when all points are allowed to freely vary. This is an 



this problem is beyond the scope of this review and will 



appear elsewhere (Chen, et al. 2011) 



issue for hedgehogs where, as we saw in Section III drag- 
ging a point defect around a disclination line required us 
to carefully keep track of the base point. Although at the 
end of the process the base point returned to its original 
value in EP^, the journey was not uneventful and re- 
sulted in the reversal of the sign of the hedgehog charge. 
Thus while d and —d are distinct when a base point is 
held fixed, they may be freely converted into each other 
if all points on the measuring circuit are allowed to vary. 
Without the base point the group structure is lost and 
since, as we have argued, it is more natural to do without 
the base point, the combination of defects in condensed 
matter does not normally follow the laws of group com- 
position^^; cases such as the XY or Heisenberg models 
are exceptions to this norm. Of course the group compo- 
sition can be regained, but only at the expense of keeping 
track of paths that the defects move along - the tethers 
of Section [IVl 

Throughout this colloquium we have avoided dis- 
cussing any details of how the defects themselves ap- 
pear and disappear. Indeed, when the authors have tried 
to observe liquid crystal textures under the microscope, 
we have noted, upon cooling, the sudden appearance of 
many defects followed by a coarsening of defects too rapid 
for us to nicely photograph^^. The coarsening is a situa- 
tion in which the defects become dynamical objects. Not 
only do they move, they coalesce into defect/anti-defect 
pairs and disappear. Though the addition rules afforded 
by the general homotopy group structure are preserved, 
the sample itself changes. Since defects are places in 
which the local order is ill-defined, we can count the num- 
ber of defects as well as measuring their charge. When 
defects annihilate the defect number changes! To account 
for this we have to perform "surgery" where we remove 
an arbitrarily small region of the sample along an inci- 
sion and replace it with a new region that matches the 
texture smoothly along the cut. A proper formulation of 



This is even more vivid in systems with broken translational 



symmetry such as smectic hquid crystals I Chen, et al. 2009 i 



^® This may reflect more the experimental abilities ot the authors 
than the true nature of liquid crystals. 



We have intended to bring to light the subtleties of the 
theory of topological defects in experimentally realizable 
systems. With the advent of recent experimental work 
on colloidal inclusions in liquid crystals, these issues are 
no longer purely academic. We hope we have given the 
reader the tools to go forth, combine, and multiply by 
themselves. 
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